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ABSTRACT 

A Suzaku observation of the nucleus of the radio-loud AGN Centaurus A in 2005 has yielded a broad- 
band spectrum spanning 0.3 to 250 keV. The hard X-rays are fit by two power laws, absorbed by columns 
of 1.5 and 7 x 10 23 cm -2 . The dual power-laws are consistent with previous suggestions that the power- 
law components are X-ray emission from the sub-pc VLBI jet and from Bondi accretion at the core, 
or are consistent with a partial covering interpretation. The soft band is dominated by thermal emis- 
sion from the diffuse plasma and is fit well by a two-temperature VAPEC model, plus a third power-law 
component to account for scattered nuclear emission, kpc-scale jet emission, and emission from X-ray 
Binaries and other point sources. Narrow fluorescent emission lines from Fe, Si, S, Ar, Ca and Ni are 
detected. The width of the Fe Ka line yields a 200 light-day lower limit on the distance from the black 
hole to the line-emitting gas. K-shell absorption edges due to Fe, Ca, and S are detected. Elemental 
abundances are constrained via the fluorescent lines strengths, absorption edge depths and the diffuse 
plasma emission lines. The high metallicity ([Fe/H]=+0.1) of the circumnuclear material compared to 
that in the metal-poor outer halo suggests that the accreting material could not have originated in the 
outer halo unless enrichment by local star formation has occurred. Relative abundances are consistent 
with enrichment from Type II and la supernovae. 

Subject headings: galaxies: active — X-rays: galaxies — galaxies: individual (NGC 5128) 


1. introduction 

The radio-loud active galactic nucleus (AGN) Centau- 

rus A (NGC 5128) is one of the most extensively-studied 

AGNs at all wavebands, thanks to its proximity (distance 

of 3.8 ± 0.4 Mpc, Rejkube 2004; 1' = 1 kpc) and bright- 

ness (2-10 keV flux typically ~2 xlO -10 erg cm -2 s -1 ). 

Cen A features the nearest AGN jet and is considered a 

prototypical FR I radio galaxy. Its radio structure con- 

sists of extended kpc-scale outer, middle, and inner radio 

lobes, fed by a one-sided kpc-scale jet. The jet is at an 

angle of ~60° to the line of sight, leading to its descrip- 
tion as a ’misaligned’ BL Lac object (e.g., Bailey et al. 


1986). VLBI observations have also revealed a sub-pc jet 
(see Israel 1998 for a review). 

The host galaxy of Cen A is a giant elliptical with a pro- 
nounced optical dust lane that obscures the inner few kpc 
and is believed to be an edge-on disk structure (Quillen et 
al. 1992). The inner dust disk, with associated young stars, 
H II regions, and other ionized gas, is the site of vigorous 
star formation, e.g., as observed along the edges of the in- 
ner disk (e.g., Ebneter & Balick 1983, Dufour et al. 1979). 
The dust disk is interpreted as the result of a merger be- 
tween a large elliptical and a gas-rich disk galaxy at least 
a Gyr ago (e.g., Israel 1998). Infrared studies are neces- 
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sary to reveal the compact, pc-scale core (e.g., Karovska 
et al. 2003). NIR studies (e.g, Schreier et al. 1998) have 
revealed a circumnuclear disk tens-hundreds of pc wide, 
surrounding the central black hole. Stellar kinematic stud- 
ies indicate a black hole mass near 2xl0 8 Mq (Silge et al. 
2005) 

Previous studies of the nuclear X-ray emission from ~4 
keV to several hundred keV (e.g., Rothschild et al. 1999) 
have established the presence of a non-thermal, power- 
law continuum whose origin is not certain. It could be a 
signature of accretion, or it could be associated with jet 
emission processes, e.g., synchrotron or inverse Compton 
emission from the sub-pc VLBI jet. Below ~4 keV, the 
X-ray continuum emission undergoes moderately heavily 
absorption, though the nature of the X-ray obscuring ma- 
terial is not certain, i.e., it could have a sky-covering frac- 
tion of 1 as seen from the central X-ray source or it could 
be in a form of a dusty torus. Rothschild et al. (2006) have 
noted ~50% variations in the absorbing column over ~20 
years, and suggested that the nucleus is seen through the 
edges of a warped, rotating disk. 

Previous hard X-ray spectral fits to the nucleus of Cen 
A center on one or more power-laws, with varying degrees 
of absorption. For instance, using ROSAT and ASCA, 
Turner et al. (1997) suggested that one plausible descrip- 
tion of the nuclear emission is a partial-covering model, 
wherein the nucleus is seen through three layers of absorp- 
tion: 40%, 59%, and 1% of the nuclear emission is obscured 
by columns of 4, 1 and 0.01 x 10 23 cm -2 , respectively. 
Using Chandra- HETGS and XMM-Newton observations, 
Evans et al. (2004) suggested a different interpretation. 
They fit the hard X-ray continuum spectrum by a heavily- 
absorbed power-law, thought to be associated with Bondi 
accretion at the core, plus a relatively less-absorbed power- 
law thought to be associated with the sub-pc VLBI jet, 
though they did not rule out the partial covering interpre- 
tation. 

There is a prominent Fe Ka emission line at 6.4 keV, 
first noted by Mushotzky et al. (1978). It is narrow (e.g., 
FWHM = 2200±900 km s" 1 , Evans et al. 2004), and its 
flux is historically nearly constant despite long-term con- 
tinuum variations spanning a factor of ~5 over 20 years 
(Rothschild et al. 1999), implying an origin for the line 
which is distant from the origin of the variable continuum. 

The soft X-ray emission in Cen A is a combination of 
narrow emission lines plus a continuum of some form, and 
has generally been interpreted as thermal emission, likely 
associated with the diffuse gas that extends to ~6 kpc 
around the nucleus. Turner et al. (1997) modeled the soft 
X-rays with a 0.6 keV thermal plasma component, a 5 keV 
component associated with emission from a population of 
X-ray Binaries, and a power-law jet component. 

In this paper, we report on an observation of Cen A 
made with the Suzaku observatory in August 2005. The 
combination of the X-ray Imaging Spectrometer (XIS) 
CCD and the Hard X-ray Detector (HXD) instruments 
have yielded a broadband spectrum covering 0.3 to 250 
keV, allowing us to deconvolve the various broadband 
emitting and absorbing components in this object. Fur- 
thermore, the exceptional response of the XIS CCD and 
high signal-to-noise ratio of this observation allows us to 
study narrow emission lines in great detail. This observa- 


tion has also yielded the highest quality soft X-ray spec- 
trum of Cen A obtained so far. §2 gives a brief overview of 
the Suzaku observatory, and describes the observation and 
data reduction. §3 describes the spectral fits. The results 
are discussed in §4, and a brief summary is given in §5. 

2. OBSERVATIONS AND DATA REDUCTION 

The nucleus of Cen A was observed by Suzaku from 2005 
August 19 at 03:30 UT until August 20 at 09:50 UT, and 
in fact was the first light target for the HXD. Suzaku was 
launched 2005 July 10 into a low-Earth orbit. It has four 
X-ray telescopes (XRTs; Serlemitsos et al. 2007), each with 
a spatial resolution of 2' (HPD). The XRTs focus X-rays 
onto four X-ray Imaging Spectrometer (XIS; Koyama et 
al. 2007) CCDs, which are sensitive to 0.2-12 keV X-rays 
on a 18' by 18' field of view, contain 1024 by 1024 pixel 
rows each, and feature an energy resolution of ~140 eV at 
6 keV. Three CCDs (XIS-0, 2 and 3) are front-illuminated 
(FI), the fourth (XIS-1) is back-illuminated (BI) and fea- 
tures an enhanced soft X-ray response. The XRT/XIS 
combination yields effective areas per detector of roughly 
330 cm 2 (FI) or 370 cm 2 (BI) at 1.5 keV, and 160 cm 2 (FI) 
or 110 cm 2 (BI) at 8 keV. Each XIS is equipped with two 
55 Fe calibration sources which produce fluorescent Mn K a 
and K f3 lines and are located at the CCD corners. Suzaku 
also features a non-imaging, collimated Hard X-ray Detec- 
tor (HXD; Takahashi et al. 2007); its two detectors, PIN 
and GSO, combine to yield sensitivity from ~10 to ~700 
keV. Further details of the Suzaku observatory are given 
in Mitsuda et al. (2007). 

2.1. XIS reduction 

The XIS data used in this paper were version 0.7 of 
the screened data (Fujimoto et al. 2007) provided by the 
Suzaku team. The screening is based on the following 
criteria: grade 0, 2, 3, 4, and 6 events were used, the 
CLEANSIS script was used to remove hot or flickering pix- 
els, data collected within 256 s of passage through the 
South Atlantic Anomaly (SAA) were discarded, and data 
were selected to be 5° in elevation above the Earth rim 
(20° above the day-Earth rim). The XIS-FI CCDs were in 
2x2, 3x3 and 5x5 editmodes, for a net exposure time after 
screening of 71.5 (XIS-0), 70.3 (XIS-2) and 68.6 (XIS-3) 
ks. XIS-1 was in 3x3 and 5x5 modes, for a net exposure 
of 68.6 ks. 

The source was observed at the nominal center position 
of the XIS. To reduce the risk of photon-pileup, the XIS 
was kept in “1 /4 Window Mode” for this observation: each 
CCD recorded 256 by 1024 pixel rows, with a readout time 
of 2 s, as opposed to 8 s when the full 1024 by 1024 win- 
dow is used (“Normal Mode”). For each XIS, we extracted 
a 2' radius centered on the source. Spectra were binned 
to a minimum of 50 counts bin -1 to allow use of the x 2 
statistic. 

Because the soft diffuse emission in Cen A covers nearly 
the entire read-out area of the CCD, we used observations 
of the North Ecliptic Pole (NEP) to extract a background 
spectrum. Suzaku observed the NEP on 2005 September 
2 to 4. There was a brief, sharp flare during the NEP ob- 
servation, consistent with evidence for solar wind charge 
exchange emission (see Fujimoto et al. 2007). Removing 
that flare and applying the same screening criteria as above 
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yielded a net exposure time of 93.9 (XIS-0), 108.5 (XIS-1), 
93.3 (XIS-2), 96.6 (XIS-3) ks. We extracted data over a 2' 
circle. The NEP is known to contain soft X-ray emission 
lines (see e.g., Fujimoto et al. 2007), but they are much 
fainter than those of Cen A (see §3). The O VII and O 
VIII lines in Cen A are a factor of roughly 10 (20) higher 
compared to NEP as seen in the XIS-BI (XIS-FIs). Lines 
near 0.8-0. 9 keV such as Fe L XVII and Ne IX in Cen A 
are typically ~60 (100) times brighter in Cen A compared 
to the NEP as seen in the XIS-BI (XIS-FIs). 

Average 0.5-2.0 keV count rates for the NEP observa- 
tion were 0.004, 0.009, 0.004 and 0.003 count s -1 for XIS-0, 
1, 2 and 3, respectively; average 2-10 keV count rates were 
0.005, 0.012, 0.005 and 0.004 count s'" 1 , respectively. 

Response matrices and ancillary response files (ARFs) 
were generated for each XIS independently using XISSIM- 
RMFGEN and XISSIMARFGEN version 2006-10-26 (Ishisaki 
et al. 2007). The ARF generator takes into account the 
level of hydrocarbon contamination on the optical block- 
ing filter. However, the Cen A observation occurred very 
early in the Suzaku mission, and the level of contamina- 
tion was quite low: we estimate a carbon column density 
of only ~0.1, 0.2, 0.3 and 0.5xl0 18 cm -2 for XIS-0, 1, 2 
and 3, respectively. 

A consequence of observing in 1 /4 Window Mode is that 
there are no 55 Fe calibration source lines on the CCD. 
This leads to a higher systematic uncertainty in the en- 
ergy scale compared to normal mode (where it is a few 
tenths of a percent at most). To estimate the instrument 
resolution, we used the calibration source lines from a 77 
ks observation of MCG-6-30-15, observed in Normal Mode 
immediately prior to Cen A. We fit the calibration source 
spectra with three Gaussians. Two Gaussians were for the 
Mn Ka doublet (expected energies 5.899 keV and 5.888 
keV), with energy centroids fixed to be 11 eV apart, and 
the higher energy line flux set to twice that of the lower 
energy one. The third Gaussian was used to model the K/3 
line, expected at 6.404 keV. We found the average of all 
the calibration line widths a to be 9^2 eV. 

The positional accuracy of Suzaku is ~1' at present, as 
the spacecraft has been known to exhibit small attitude 
variations (“wobble”). We checked to make sure that the 
effect of these attitude variations was not significant, given 
that the Cen A observation was in 1/4 Window Mode. We 
generated light curves by extracting over radii of 1.5' and 
1.0' and looked for variations on the orbital timescale, but 
found nothing significant. 

Average 0. 5-2.0 keV net source count rates were 0.24, 
0.35, 0.25, and 0.25 count s _1 for XIS-0, 1, 2, and 3, re- 
spectively. Average 2-10 net source count rates were 4.7, 
4.4, 4.5, and 4.5 count s -1 for XIS-0, 1, 2, and 3, re- 
spectively. Figure 1 shows 2-10 keV light curves for each 
XIS. During the observation, the source flux increased by 
only 10%, with fractional variability amplitudes F var (see 
Vaughan et al. 2003 for definition) of 3.7 ±0.2 %, 3.4 ±0.2 
%, 3.1 ± 0.2 %, and 3.3 ± 0.2 % for XIS-0, 1, 2, and 3, re- 
spectively. The 0. 5-2.0 keV flux (determined from model 
fits; see §3) was 2.33xl0~ 12 erg cm -2 s" 1 . The 2-10 keV 
flux was 2.12xlO“ 10 erg cm -2 s -1 . 


2.2. HXD reduction 

The PIN and GSO source spectra were extracted from 
cleaned version 1.2 (pre-1.2-rl) HXD event files provided 
by the HXD instrument team. We first discuss the PIN ex- 
traction. Data were selected according to the following cri- 
teria: at least 500 s since SAA passage, COR > 8 GV, and 
day- and night-Earth elevation angles each >5°. Instru- 
mental (non-X-ray) background spectra for the PIN were 
provided by the HXD Team (“Background D” model) gen- 
erated from a time-dependent model. 12 The model utilized 
the count rate of upper discriminators as the measure of 
cosmic ray flux that passed through the silicon PIN diode 
and yielded background spectra based on a database of non 
X-ray background observations with the PIN (Fukazawa 
et al. 2007). The systematic uncertainty of the PIN is 
expected to be <5%. Both the source and background 
spectra were generated with identical good time intervals, 
and the exposures were corrected for instrument dead time 
(a ~5% effect). This yielded a good time exposure of 60.8 
ks. Data <12 keV were discarded due to noise contam- 
ination near the lower threshold of the PIN diode. Data 
above 76 keV were also discarded: the gain above an in- 
ternal Bi Ka calibration line at 76 keV is not well-defined, 
though there is GSO data covering these energies. Fur- 
ther details of the HXD in-orbit performance are given 
in Kokubun et al. (2007). To model the contribution to 
the total background from the Cosmic X-ray Background 
(CXB), the spectrum of the form 9.0xl0“ 9 (E/3keV)“°‘ 29 
exp(-E/40keV) erg cm" 2 s _1 sr” 1 keV -1 (Gruber et al. 
1999) was used; the contribution in the 12-76 keV band 
was 1.1 xlO -11 erg cm”" 2 s _1 . 

The spectrum was binned to a minimum of 
400 count bin"" 1 . We used the response file 
aeJixd_pinxinom_20060814.rsp. The 12-76 keV net source 
flux and count rate were 7.3xlO“ 10 erg cm -2 s _1 and 
1.21 count s'" 1 . The total (X-ray plus particle) back- 
ground 12-76 keV flux and count rate were 6.9xlO“ 10 
erg cm -2 s” 1 and 0.59 count s -1 . The or bitally-binned 
light curve is shown in Figure 1. The light curve gener- 
ally increases in flux by ~ 10-20% over the observation, 
matching the general trend shown by the XIS light curve, 
though shorter-timescale differences between the XIS and 
PIN light curves might be attributed to uncertainty as- 
sociated with the PIN background. Figure 2 shows the 
net source, background, and total (source ± background) 
spectra. The source spectrum is always at least 30% of 
the total up to ~50 keV. 

The HXD-GSO data were reduced in a manner simi- 
lar to the HXD-PIN data. A key contributor to the GSO 
particle background is activation lines, e.g., delayed emis- 
sion from radioactive isotopes induced inside the detector 
due to interaction with SAA particles. To minimize such 
emission, data taken within 6000 s of SAA passage were 
discarded; the background is the most reliable and the 
least variable during non-SAA orbits. This yielded a good 
time exposure of 17.1 ks. The ’Background D’ model files 
provided by the HXD Team were used (Fukazawa et al. 
2007). The CXB is expected to contribute insignificantly 
to the total GSO background and was ignored. Source 


12 Other recently-published Suzaku results have made use of the “Background A” model. However, observations occurring before September 
2, 2005 had a hit-pattern width set to a shorter value compared to observations occurring after this date. Because applying the “Background 
A” model would yield an underestimate of the true non-X-ray PIN background, we instead use the “Background D” model. 
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and background spectra were both binned with GRPPHA 
following recommendations from the HXD Team 13 . 

The source was detected out to 250 keV; the GSO is 
sensitive down to roughly 45 keV. The 45-250 keV net 
source flux and count rate were 7.2 xlO -10 erg cm -2 s _1 
and 0.83 count s -1 . The 45-250 keV background flux and 
count rate were l.OxlO -8 erg cm -2 s -1 and 9.1 count 
s -1 . Figure 1 shows the GSO light curve. Figure 2 
shows the net source, background, and total (source 4- 
background) spectra. The source spectrum is always at 
least 5% of the total below 200 keV. The response file 
aeJixd_gsoxinom_20060321.rmf was used. 

Below 100 keV, the GSO field of view is 34' x 34' FWHM, 
the same as the PIN. Above 100 keV, the field of view in- 
creases with photon energy to a maximum of 4.5° x4. 5° 
FWHM and there is the possibility of source confusion. 
The nearest possible contaminating source is the blazar 
MS 1312.1-4221 = GRO J1312-42, located ~2° west of 
Cen A. According to Steinle et al. (1998), the flux of 
MS 1312.1-4221 in the Compton Gamma- Ray Obervatory 
( CGRO)-OSSE band (50-4000 keV) is likely less than a 
tenth of Cen A’s flux. MS 1312.1-4221 is much fainter at 
other X-ray bands as well: the 1 keV intensity is usually 
only 1-10% of the historically-observed Cen A intensities 
(Kinzer et al. 1995). We henceforth assume that the con- 
tamination from MS 1312.1-4221 is negligible in the GSO 
band. 

3. SPECTRAL ANALYSIS 

Given the 2' FWHM HPD of the XRT, the spectrum 
will contain blended contributions from multiple sources. 
Previous studies have shown the hard X-rays 14 to be dom- 
inated by two power-laws; both hard X-ray sources are 
constrained to lie within a few tenths of an arcsecond 
(~5 pc) of the position of the unresolved radio core by 
Chandrar ACIS (Kraft et al. 2000). The soft X-rays con- 
sist of thermal emission from the diffuse plasma and dif- 
fuse emission from the kpc-scale jet, with contributions 
from a few of the innermost knots in the kpc-scale jet and 
other point sources resolved by Chandra- ACIS (Kraft et 
al. 2000; Kataoka et al. 2006). Each knot, however, is ex- 
tremely faint in the soft X-rays (< 6 x 10" 13 erg cm" 2 
s -1 ; Feigelson et al. 1981). Evans et al. (2004) noted that 
the 1 keV flux density of the power-law emission associ- 
ated with the sub-pc jet was a factor of ~15 brighter than 
that associated with the kpc-scale jet. 

The four XIS spectra were fit separately. 0.4-11.5 keV 
FI data and 0.3-11.5 keV BI data were included. We found 
it reasonable to keep the photon indices tied for all four 
XlSes while allowing the relative normalizations for XIS- 
0 and XIS-1 to each vary relative to that for XIS- 2 and 
3, which were tied together (XIS-3/XIS-2 fixed at 1.00). 
XIS-O/XIS-2 was typically 1.06; XIS-l/XIS-2 was typically 
0.96 (in our best-fit model, XIS-O/XIS-2 = 1.065dh0.004 
and XIS-l/XIS-2 = 0.951±0.003). We ignored 1.80-1.87 
keV in the XIS FI spectra and 1.75-1.87 keV in the XIS 
BI spectrum due to uncertainties in calibration associated 
with the instrumental Si K edge. This means we cannot 


directly study the intrinsic neutral Si K edge in the ab- 
sorbing material or any highly-ionized Si K lines expected 
due to the thermal plasma emission. The PIN /XIS-2 and 
GSO/XIS-2 normalizations were left free, but were usu- 
ally close to 1.09 and 0. 9-1.0, respectively, though results 
were not strongly dependent on these factors (in our best- 
fit model, PIN /XIS-2 = 1.055±0.022 and GSO/XIS-2 = 
0.98±0.06). Preliminary results on the relative GSO/PIN 
normalization using the Crab suggest GSO/PIN = 0.84 
(Yamasaki et al., in prep). Fixing GSO/PIN at 0.84 in 
the Cen A fits usually resulted in x 2 increasing by <5; 
the photon index did not change significantly. We left the 
GSO /PIN normalization free to attain the lowest x 2 possi- 
ble. All errors on one interesting parameter correspond to 
Ax 2 = 2.71 (with all relative normalizations except XIS- 
3/XIS-2 left free). A redshift of 0.001825 (Graham 1978) 
was used. Galactic absorption of 8.6 xlO 20 cm -2 was in- 
cluded. The abundances of Lodders (2003) were used. Ab- 
solute abundances relative to solar values are denoted by 
Z Fe ; a value of 1 denotes solar abundance. 

As shown below, the best model fit includes five over- 
lapping broadband components, and changes in how one 
is fit can potentially cause changes in other components. 
Our strategy in finding the best-fit broadband model was 
to fit the hard band first, then extend it to the soft band. 
Relatively more narrow features (lines and edges) were ex- 
plored later. 

3.1. Broadband fit 

We first fit the >3 keV emission with an absorbed power- 
law (henceforth denoted PL1) using ZVPHABS(cutoffpl) 
(Model 1). The power-law cutoff was kept fixed at 400 
keV; abundances in the absorber were initially kept fixed 
at solar. Residuals are plotted in Figure 3. The slight 
hard excesses above ~9 keV in the XIS spectra are asso- 
ciated with calibration uncertainties known at the time of 
this writing. Residuals in the GSO spectrum near 150-200 
keV are likely associated with uncertainties in estimating 
the strengths of activation lines in the GSO background 
(lines at 153 and 196 keV are due to 153 Gd and 151m Eu, 
respectively); see Kokubun et al. (2007) for further details. 
GSO residuals near 65 keV are likely associated with cal- 
ibration uncertainties present at the time of this writing. 
Large residuals near the expected Fe Ka line were appar- 
ent. We then added Fe Ka and K/3 lines; the K/3 energy 
was fixed at 7.056 keV, its intensity was fixed at 0.13 times 
that of Ka, and the widths of the two lines were tied, x 2 
dropped from 1.50 to 1.07, but the residuals in the XIS 
band suggested the fit could be improved if we added a 
second, fainter absorbed power-law (PL2; Model 2). The 
photon index of the second power-law, T 2 , was kept tied to 
that of PL1, T i, for simplicity, but the absorbing column 
densities Nm and A^h 2 were not. x 2 dropped by over 140. 
PL1, with a 1 keV normalization ~4 times that of PL2, was 
absorbed by a column JVhi near l.OxlO 23 cm -2 ; Ah 2 was 
near 5.0xl0 23 cm -2 . The photon index was 1.773+0. oil- 
Residuals to Model 2 are plotted in Figure 3. 


13 group 0 24 25 25 26 2 27 28 2 29 31 3 32 35 4 36 38 3; group 39 42 4 43 46 4 47 51 5 52 56 5 57 62 6; group 63 68 6 69 75 7 76 83 8 84 91 8 
92 100 9; group 101 110 10 111 121 11 122 134 13 135 147 13 148 162 15; group 163 178 16 179 196 18 197 216 20 217 238 22 239 262 24; group 
263 288 26 289 317 29 318 349 32 350 384 35 385 422 38; group 423 465 43 466 511 46 

14 We henceforth denote “hard X-ray” emission as any intrinsic emission component above 2-3 keV, as opposed to describing emission detected 
solely with the HXD. 
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We then included the soft band data, down to 0.3 keV. 
Several soft X-ray emission lines are prominent, and are 
especially well constrained with XIS-1, as illustrated in 
Figure 4. Simply adding a third power-law (Model 3) to 
the soft X-rays thus did not yield an acceptable broadband 
fit. Table 1 shows results for the lines when each is fit with 
a Gaussian; we identify them as originating in O VII, O 

VIII, a blend of Fe XVII 3s-2p and O VII RRC, a blend 
of Fe XVII 3d-2p ('PJ and Fe XVII 3d-2p ( 3 Di), Ne 

IX, Ne X and Mg XI. We can use the lines as a check of 
the absolute energy scale of the XIS. Assuming that the 
He-like lines are each dominated by the forbidden line, as 
expected for a thermal plasma, the measured energy cen- 
troids, including Ne X are within <10 eV of the expected 
energies (within <4 eV excluding Ne X). Additionally, the 
Fe K a centroid is 6±3 eV from its expected energy. 

Attempting to model the soft band emission with 
a single-temperature VAPEC component, with all abun- 
dances fixed at solar, also did not yield an acceptable fit. 
In Model 4, we tried two VAPEC components, henceforth 
denoted VAPEC 1 and VAPEC2. The best-fit model had 
temperatures of k^T ~ 0.2 keV, to model the He- and H- 
like O lines, and k^T ~ 0.6 keV, to model the Fe L and 
He- and H-like Ne lines. However, x 2 was still high? about 
1.10, and there were still large residuals <0.7 keV. 

It is plausible that some optically-thin circumnuclear 
material may scatter some of the hard X-ray continuum; a 
scattered power-law component was included in the best- 
fit model of Turner et al. (2005), for instance. In addition 
to the two VAPEC components, we added a third power- 
law (PL3), with photon index T 3 tied to T i, assuming that 
PL3 undergoes no absorption in excess of the Galactic col- 
umn. x 2 /dof was high: 10158/9161. Untying T 3 improved 
the fit considerably: x 2 1 dof fell to 9806.8/9158; T 3 was 
1.16±g;Jf (Model 5). 

However, it is plausible that PL3 does undergo some ab- 
sorption in the host galaxy, so we added an ZVPHABS com- 
ponent to PL3, and refit with T 3 tied to IT. x 2 /dof fell 
to 9798.9/9160. The column density of the third absorber, 
7 Vh 3> was 0.13±0.06 xlO 22 cm -2 . The 1 keV normaliza- 
tion of PL3 was 0.8% of the sum of the normalizations of 
PL1 and PL2. Untying T 3 improved the fit considerably: 
X 2 /dof fell to 9745.1/9159, T 3 was 1.28;to.i2> an< ^ ^H3 was 
<0.03 xlO 22 cm -2 (Model 6). All of the improvement in 
the fit when T 3 was thawed came at energies below 2 keV. 
PL3 likely represents a blend of scattered emission (likely 
< 0.8% of the total nuclear hard X-ray continuum) plus 
X-ray emission associated with the kpc-scale jet, with per- 
haps some contribution from XRBs or ULXs in the inner 2 
kpc of the host galaxy. For the moment, we adopt Model 
6 as our new baseline broadband model. Best-fit parame- 
ters for Models 5 and 6 are listed in Table 2. Data/model 
residuals for Models 3, 4, and 6 are shown in Figure 5. 

3.1.1. Hard X-ray continuum emission 

Evans et al. (2004) fit the hard X-ray Chandra and 
X MM- Newton data with a dual power-law model, wherein 
T i and T 2 differed by about 0.3. We untied T 2 from T\. 
However, compared to Model 6, x 2 only dropped by 1.4, 
and the best-fit values of Ti and IT were consistent; the 
uncertainty on V 2 was ±0.24. In subsequent fits, we will 
continue to fix T 2 =T\. 


There has been no strong evidence from past observa- 
tions for a strong Compton reflection hump in Cen A. We 
added a Compton reflection component to Model 6 using 
PEXRAV, keeping the input normalization equal to that 
of PL 1, the inclination fixed at 30°, the high-energy cutoff 
fixed at 400 keV, assuming solar abundances, and leaving 
the PIN/XIS-2 normalization free. The best fit-model pre- 
ferred a value of the reflection fraction R (defined as fi/27r, 
where is the solid angle subtended by the reflector) of 
0, with an upper limit of 0.05, consistent with previous 
suggestions that the Fe Ka line originates in Compton- 
thin material. This limit is identical to that obtained by 
Benlloch et al. (2001) using Rossi X-ray Timing Explorer 
( RXTE) data. Contour plots of R versus T i and versus the 
PIN/XIS-2 relative instrument normalization are shown in 
Figure 6. 

Using CGRO- OSSE, Kinzer et al. (1995) found evidence 
for a high-energy cutoff which appeared at 300 keV to 
~700 keV in the highest- to lowest-flux states, respec- 
tively. Keeping the cutoff energies for PLl and PL2 in 
Model 6 equal to each other, we found a lower limit of 
400 keV. Future improvements to background modeling of 
the GSO may eventually allow us to establish more strict 
constraints. 

3.2. Hard X-ray emission lines and absorption edges 

In addition to that from Fe, fluorescent Ka lines from 
Si and S have been claimed by Sugizaki et al. (1997) and 
Evans et al. (2004). We added five more Gaussians, each 
with width tied to that of the Fe Ka line; it was signifi- 
cant at >90% confidence according to the F-test to add 
each line. The centroid energies are consistent with Ka: 
emission from Si, S, Ca, Ar, and Ni; detection of the lat- 
ter three lines are reported here for the first time. The 
decrease in x 2 for each line, including the Fe Ka and K/3 
lines, along with energy centroids, fluxes, and observed 
equivalent widths FWobs? are listed in Table 3. In the case 
of Si, we caution that, as data near 1.8 keV have been ig- 
nored due to calibration uncertainties, values of flux and 
equivalent width may be uncertain (e.g., some fraction of 
the large value of Ax 2 may be due to calibration uncer- 
tainties). The detection of the Ni line is robust: there 
does exist a Ni Ka line in the NEP background spectrum, 
but it is 200 times fainter than the source spectrum line. 
Data/model residuals to models with all six Ka lines and 
the Fe K/3 line removed are shown in Figures 7-9. 

The best-fit Fe Ka line width is 14 eV, but formally the 
value is only an upper limit, 24 eV. After subtracting the 
55 Fe calibration line width of 9j^ eV in quadrature, we 
infer an intrinsic line width of < 23 eV, which corresponds 
to a FWHM velocity of <2500 km s -1 . 

We next explored the abundances of the hard X-ray ab- 
sorbers (Model 7). The K-shell edges for Fe, Ca, and S are 
detected at >99.999% confidence: setting each ZVPHABS 
abundance to 0. The detection of the Ca K edge is likely 
robust: the response of the Fe Ka line includes a Si es- 
cape feature near 4.5 keV, but it is a factor of more than 
50 fainter than the continuum and likely does not influence 
the Ca edge. The intrinsically weak Ar K and Ni K edges 
are not detected at high confidence; in the latter case, the 
decrease in the effective area of the XIS above 8 keV is also 
a factor. We re-fit the data, allowing Zp e , Zca. and Zs to 
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vary, and keeping the abundances of all three absorbers 
tied to each other for simplicity. We also kept all other 
abundances fixed at solar. Best-fit model parameters for 
Model 7 are listed in Table 2; the edges are also displayed 
in Figures 7-9. The abundance results are listed in Table 
4; note that the Fe abundance is inconsistent with solar at 
3cr confidence. Given the column density of the absorbing 
material and the 2-3 keV rollover, most of the rollover is 
the result of absorption by K-shell O (with some contribu- 
tions to the total opacity from K-shell Ne and L-shell Fe). 
We are actually measuring abundances relative primarily 
to O, not to H. The abundances listed relative to H in 
Table 4 implicitly assume Zo = 1.0. 

To find the equivalent optical depth for each edge, we 
re-fit the data with Zp e , Zc a and Zs set to zero in the 
ZVPHABS components, but with edges at the appropriate 
energies. The results, listed in Table 4, show that the en- 
ergy of each edge is consistent with absorption by neutral 
atoms. 

The Fe K/3 /Ka intensity ratio can yield insight into the 
ionization state of the line-emitting gas. However, in the 
present spectrum, the K/3 line is somewhat blended with 
the Fe K edge. However, we note that the Fe K edge energy 
and depth are robust to the properties of the K/3 line, and 
the energy can be used to constrain the ionization state of 
the material, as discussed further in §4. 

We searched for any possible Compton shoulder to the 
Fe Ka line by adding a Gaussian near 6.24 keV, but there 
was no improvement in the fit. We find an upper limit 
of 1.7xl0“ 5 ph cm -2 s -1 to Compton shoulder emission 
(EW < 10 eV, determined relative to locally-fit contin- 
uum). 

Grandi et al. (2003) claimed emission due to ionized Fe 
near 6.8 keV using BeppoSAX. We added a narrow Gaus- 
sian at the rest energies of Fe XXV and Fe XXVI, with 
no change in the fit statistic in either case; upper limits to 
emission were 3 and 1 eV, respectively. 

3.3. Abundances of the thermal plasma 

We next explored the abundances Z of the diffuse, ther- 
mal plasma relative to solar values. Abundances for the 
VAPEC1 and VAPEC2 components were kept tied to each 
other. We thawed one element at a time from solar, while 
assuming solar abundances for all other elements. Abun- 
dances relative to H can be measured accurately if the 
continuum is well-determined and the thermal emission is 
not too heavily dominated by lines. Given the presence of 
PL3, one must keep in mind that there could be systematic 
effects associated with some derived abundance values. 

We did not find any strong evidence for non-solar abun- 
dances for C, N or O. The best chance to constrain Zq 
in our model is via the C VI Ly a and Ly /3 lines at 
0.368 and 0.436 keV, respectively. However, there are 
very few counts at these energies, and the PL3 component 
dominates at these energies. Similarly, Z n is constrained 
mainly via the N VII Lya line, which is not robustly de- 
tected. Allowing Zo to vary resulted in a best-fit value 
of 0.95±0.35 solar, with x 2 decreasing by only 1.7. We 
assume Zq, Z n and Zo are 1 for the rest of the paper. 

We did not test the Si abundance: Si XIII emission can- 
not be constrained, as we have ignored data near the Si K 
instrument edge. The best-fit model does predict a modest 


Si XIV emission line at 2.05 keV, but at that energy, the 
PL3 component and the absorbed hard X-ray power-laws 
dominate. 

In addition to Fe XVII emission lines, Fe is also respon- 
sible for a large bump in the continuum emission from 
about 0.7 to 1.0 keV. Adjusting the Fe abundance causes 
the PL3 component to vary in normalization to compen- 
sate, creating systematic uncertainty to Zp e derived from 
L-shell emission. We henceforth rely on the abundance 
determined using the K edge depth (§3.3). 

We did, however, find significant improvements to the 
fit when the Ne or Mg abundances were thawed. In the 
case of Ne, the best-fit abundance was 2. 71^2 s °lar; X 2 de- 
creased by 71 (F-test probability 2xl0“ 15 ). For Mg, the 
best-fit abundance was 1.6±0.3 solar, with x 2 decreasing 
by 19 (F-test probability 2xl0” 5 ). 

We re-fit the data with the Ne and Mg abundances in 
VAPEC1 and VAPEC2 thawed (Model 8). We left the 
VAPEC abundances for C, N, O, Si and Fe fixed at solar 
(fixing Zp e at 1.17, the value derived from the Fe K edge, 
had minimal impact). In the best-fit model, Z^ e and Zug 
were 2.7±0.3, and 1.9±0.4, respectively. The tempera- 
ture of the VAPEC2 component fell slightly to 0.56io!o3 
keV. With the exception of these parameters, all other 
previously-derived parameters did not change significantly 
compared to Model 7. The best-fit parameters for Model 
8 are listed in Table 2. A contour plot of Nm versus the 
photon index and normalization of the primary power-law 
is shown in Figure 10. An unfolded model spectrum in 
plotted in Figure 11. 

4. DISCUSSION 

The broad bandpass of the XIS and HXD aboard Suzaku 
has allowed us to obtain a high-quality spectrum spanning 
nearly 3 decades in photon energy. The long exposure 
time, narrow response, and the low background make the 
XIS spectrum one of the highest quality CCD spectra of 
Cen A to date, particularly below 2 keV. 

In §4.1, we compare our results to other recent works 
to test for evidence of long-term variations in several ob- 
served parameters. The dual power-laws which dominate 
the hard X-rays, and the material which absorbs them, 
are both discussed in §4.2. That same absorbing mate- 
rial is likely the origin of the fluorescent emission lines; 
the nature of the Fe Ka line is discussed in §4.3. The 
soft X-ray emission from the thermal plasma is discussed 
in §4.4. Finally, in §4.5, element abundances are derived 
from the hard X-ray absorption edges and fluorescent emis- 
sion lines, the first time abundances in the absorbing gas 
of Cen A have been measured. In conjunction with abun- 
dances measured using the thermal plasma lines, we show 
that the abundances are consistent with enrichment from 
star formation. 

4.1. Comparison to previous observations 

We can compare our results to previous works to search 
for any variations in flux, photon index, column densities, 
and Fe Ka line flux. 

We find an absorbed 2-10 keV flux of 2.12X10" -10 erg 
cm -2 s" 1 , similar to values obtained by Rothschild et al. 
(2006) for RXTE and INTEGRAL observations of Cen A 
between 2003 and 2004. We find a 20-100 keV flux of 
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6.4x10“ 10 erg cm -2 s” 1 , within the range of 20-100 keV 
fluxes measured by Rothschild et al. (2006). We measure 
a unabsorbed 4-7 keV continuum flux of 1.5x10“ 10 erg 
cm“ 2 s _1 and a 2-10 keV absorbed luminosity of 2.9 xlO 41 
erg s” 1 , values almost identical to those found by Evans 
et al. (2004) during the 2001 XMM-Newton observation. 
These facts suggest the source has not varied in hard X-ray 
flux greatly over the last five years. 

The photon index, power-law normalization, and ab- 
sorbing column density of the primary power-law are all 
similar to values obtained by Rothschild et al. (2006). 
Evans et al. (2004), modeling XMM-Newton observations 
in 2001 and 2002 and a Chandra observation in 2002, fit 
partial-covering models (T 1 =Y 2 ) and dual power-law mod- 
els (T 2 not tied to IT). In all cases, the parameters of the 
primary (brighter) power-law and its absorber are very 
similar to those for the Suzaku observation, suggesting a 
lack of strong variations over the last five years. 

The secondary (fainter) hard X-ray power-law in our 
model has a flux of 6xl0“ n erg cm“ 2 s -1 , 3-5 times 
brighter than the component as modeled by Evans et al. 
(2004). Furthermore, we find a column density a factor of 
>10 higher. However, this is not necessarily direct proof 
that the column density has increased since 2001-2002: 
the 2-3 keV rollover in PL2 in the Suzaku fits could be in- 
fluenced by the hotter VAPEC component and PL3, while 
Moreover, Evans et al. (2004) did not include <2 keV data 
in their fits, so a direct comparison is not straightforward. 
It could be argued that our model is more suitable since 
we are using a much wider bandpass. 

The EW 0 bs of the Fe Ka line in the Suzaku observation 
is consistent with that measured from the XMM-Newton 
observations; we measure an unabsorbed Fe Ka flux of 
2.6xl0“ 12 erg cm -2 s -1 , consistent with the Chandra and 
2001 XMM-Newton observations (Evans et al. 2004). Our 
EW 0 \) S values for the Si and S lines are consistent with 
those obtained using ChandrarKETGS. The Suzaku obser- 
vation thus shows no evidence for variability of the Fe, Si 
or S Ka lines compared to the 2001-2002 observations. 

4.2. Continuum emission components 

The dual absorbed hard X-ray power-laws observed sug- 
gest two scenarios, both of which have been suggested pre- 
viously, and both of which are consistent with the Suzaku 
data. 

There could be two physically distinct sites of hard X- 
ray continuum emission, with unabsorbed 2-10 keV lumi- 
nosities of 6.7xl0 41 erg s” 1 and 1.3xl0 41 erg s _1 , each 
independently absorbed by different columns. The ab- 
sorbers are similar in column density to absorbers seen 
in some other radio galaxies (Wozniak et al. 1998). 

Evans et al. (2004) suggested that the less-absorbed 
power-law corresponds to X-ray emission from the sub-pc 
VLBI jet, while the more-absorbed power-law arises via 
Bondi accretion onto the black hole, consistent with the 
low accretion rate of Cen A (0.2% of Eddington). As far 
as jet emission is concerned, several papers in the 1980s 
(e.g., Burns et al. 1983) suggested that the X-ray emis- 
sion was synchrotron emission. However, by constructing 
the broadband SED, Chiaberge et al. (2001) concluded 
that the X-ray continuum was primarily inverse Compton 
emission, with the synchrotron component peaking in the 


far-IR and falling by 10 15 “ 16 Hz. 

The dust lane provides an optical extinction of 3-6 mag, 
which, assuming a Galactic gas/dust ratio, corresponds to 
a column Ah ~ 5-10 xlO 21 cm” 2 . The circumnuclear 
material within a few hundred pc of the black hole, how- 
ever, is believed to have a much higher column, ~10 23 
cm” 2 . Both hard X-ray absorbers are thus more likely as- 
sociated with the circumnuclear material than the optical 
dust lane. However, Evans et al. (2004) concluded that the 
brighter power-law was the more-absorbed one, contrary 
to the best-fit model for the Suzaku observation, and so it 
is not obvious which power-law component corresponds to 
emission from Bondi accretion or from the sub-pc jet. 

Alternatively, as suggested by Turner et al. (1997), a 
partial-covering model, with only one site of hard X-ray 
emission, is applicable, as Ti and F 2 are consistent. Our 
best-fit model is consistent with the idea that 84% of the 
sky as seen from the hard X-ray source is covered by a col- 
umn 1.5 xlO 23 cm” 2 , and the remaining 16% by a column 
7xl0 23 cm” 2 . 

The soft X-ray power-law, PL3, is consistent with ab- 
sorption by the optical dust lane only and is likely the sum 
of diffuse X-ray emission and knots within the kpc scale 
jet, scattered nuclear emission, and emission from other 
point sources resolved by Chandra-ACIS . The 0.5-3 keV 
luminosity of PL3 is 1 x 10 40 erg s” 1 , similar to the sum 
of the jet and knot component luminosities as reported by 
Kraft et al. (2002) and Kataoka et al. (2006). 

Finally, we comment on the absence of a strong Comp- 
ton reflection component. The low R value could of course 
indicate the absence of a large amount of Compton-thick 
material which is brightly illuminated by the primary con- 
tinuum and also efficiently radiates the reflection contin- 
uum. An alternate suggestion is that the Compton hump 
is diluted by a power-law component, e.g., from one of the 
jets. However, the large Fe line flux observed is inconsis- 
tent with this notion. In addition, Wozniak et al. (1998) 
demonstrated that spectral fitting of several radio galaxies 
which also displayed weak or non-existent Compton humps 
is inconsistent with this notion. 


4.3. The origin of the Fe Ka line 

Assuming that the line originates in gas which is in vit- 
alized orbit around the black hole, we can estimate the 
distance R from the black hole mass to the Fe Ka line- 
emitting gas. The width of the Fe Ka line corresponds 
to a FWHM velocity ufwhm of <2500 km s” 1 . Assuming 
that the velocity dispersion is related to ^fwhm as <v 2 > 
= f^pwHM (Netzer et al. 1990), assuming a black hole 
mass Mbh of 2xl0 8 M 0 , and using GMbh = ifo 2 , we 
find a lower limit to R of 6xl0 15 m = 200 light-days. The 
best-fit edge energy, 7.103±0.015 keV, is consistent with 
absorption by Fe I, though ionization stages up to Fe ~V 
cannot be ruled out (Kallman et al. 2004). 

It is plausible that the same material that absorbs the 
hard X-rays is responsible for transmitting the fluorescent 
lines. The fact that no Compton hump or 6.2 keV Comp- 
ton shoulder are seen, indicating an origin for the Fe Ka 
line in Compton-thin material, supports this notion. As 
an estimate of the Fe Ka equivalent width expected in this 
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case, we can use the following equation: 


-S'^calc — /c^/Kq-^ 


Ij 


•^Kedge 


P(E)a ph (E)N H dE 

P(Eu n e) 


( 1 ) 


This method assumes an origin in optically-thin gas which 
completely surrounds a single X-ray continuum source and 
is uniform in column density. Emission is assumed to be 
isotropic. Here, f c is the covering fraction, initially as- 
sumed to be 1.0. uj is the fluorescent yield: the value for 
Fe, 0.34, was taken from Kallman et al. (2004). /kc* is 
the fraction of photons that go into the Ka line as op- 
posed to the K/? line; this is 0.89 for Fe I. A is the number 
abundance relative to hydrogen. Initially, we assumed so- 
lar abundances, using Lodders (2003). P{E) is the spec- 
trum of the illuminating continuum at energy E; Ey me 
is the Ka emission line energy. The illuminating contin- 
uum is assumed to be the sum of both unabsorbed hard 
X-ray power-laws in Cen A, as jEWobs values were deter- 
mined relative to the total, local continuum. We note that 
the fainter power-law component cannot reproduce the en- 
tire Fe line as EW 0 \> S is too large; the primary power-law 
must be responsible for the bulk of the line. cr p h(£’) is 
the photo-ionization cross section assuming absorption by 
K-shell electrons only; all cross sections were taken from 
Veigele (1973 15 ). Nr is the column density. The two hard 
X-ray absorbing components in Model 8 have column den- 
sities 1.5 and 7 x 10 23 cm -2 , which attenuate power-laws 
with normalizations in a 5.3:1 ratio. We therefore use an 
“effective” N# column of 2.2 x 10 23 cm” 2 . 

The value of EW ca \ c is 128 eV, a factor of 1.5 higher 
than jEWobs- Assuming that Zp e = 1-17, as the Fe K 
absorption edge depth suggests, yields EW ca \ c = 153 eV. 
Possible explanations for this discrepancy include: 1) The 
initial assumption of a 100% covering fraction for the line- 
emitting gas of 1.0 may be incorrect. Assuming isotropic 
emission, f c = 0.7 (Zp e = 1.0) or 0.6 (Zp e =1.17) may 
be more appropriate. Such values would be more consis- 
tent with an origin in a disk as opposed to a shell com- 
pletely surrounding the core. 2) The illuminating X-ray 
continuum may be anisotropic. For instance, Rothschild 
et al. (2006) found that /k q was not correlated with Nu 
and suggested that the line-emitting gas may lie along the 
VLBI jet some distance from the black hole. 3) The line- 
emitting gas may be responding to an illuminating flux 
which is ~0.6 of the observed flux. This is possible, for 
instance, if the line-emitting gas is >200 light-days from 
the X-ray continuum origin, and does not lie on our line 
of sight to the continuum origin. However, it is not obvi- 
ous that this is the case, as the 2-10 keV flux measured 
by Suzaku is generally consistent with fluxes measured by 
other missions between 2001 and 2004, as noted in §4.1. 
Continued long-term monitoring of the continuum and line 
fluxes may shed light on this issue. Alternatively, the line 
of sight from the continuum origin to the line-emitting gas 
may be blocked, by a clump of gas which absorbs 40% of 
the >7 keV flux and which does not lie along our line of 
sight to the continuum. 4) It is plausible that our initial 
assumption of a uniform-density absorber is incorrect and 
the absorbers are clumpy, with higher columns lying out 
of our line of sight to the X-ray continuum source. For 


15 http://www.pa.uky.edu/~verner/photo.html 

16 http://xdb.lbl.gov 


instance, we can use Eq. 1 of Wozniak et al. (1998), which 
gives the expected Fe Ka intensity assuming a cloud with 
a column > 10 23 cm” 2 lying off the line of sight and sub- 
tending a fraction ft/47r of the sky as seen from the source, 
as a function of the K edge optical depth and the spectral 
index and normalization of the illuminating continuum. 
We find a solid angle of 0/47T ~ 0.5 yields a predicted line 
intensity compatible with what Suzaku observes. Such a 
covering fraction suggests some sort of infinite slab, such 
as an optically-thin disk. It is interesting to note that 
Evans et al. (2004) suggest that, based on the observed 
accretion rate, the Bondi accretion mode would favor a 
high efficiency flow in the form of a disk. 

4.4. Soft X-ray emission from the thermal plasma 

The Suzaku XIS has yielded the best soft X-ray spec- 
trum obtained to date for Cen A, revealing emission lines 
which likely originate in the thermal plasma that extends 
from the nucleus to a radius of ~6 kpc. At least the inner 
parts of the gas are likely heated by the AGN central en- 
gine: O’Sullivan, Ponman & Collins (2003) demonstrated 
that the gas in the central 2 kpc has a much higher tem- 
perature then the >2 to 14 kpc gas in the halo. 

To obtain a good model fit, we have used a two- 
temperature VAPEC model, with k#T near 0.2 and 0.6 
keV. However, we are likely sampling thermal emission 
from plasma spanning a range of temperatures. For ex- 
ample, the soft X-rays are likely not completely dominated 
by emission from gas with k&T above 0.6 keV, or else we 
would detect a H-like Mg line at least as intense as the He- 
like Mg line, and the He-like Ne line would be extremely 
faint. Similarly, emission from gas with k&T < 0.20 keV 
does not dominate, as the Ne X would be extremely faint. 

4.5. Constraints on abundances 

We have derived estimates of Zs , Zca and Zp e using the 
K shell edge depths, and estimates of Zq , Z^e, and Zu g 
from the thermal plasma (for the remainder of this paper, 
we assume that the abundances in the extended thermal 
plasma and in the circumnuclear absorbing material are 
identical) . 

We can also estimate abundances relative to Fe, e.g., 
Zsi/Zpe, from the Ka fluorescent lines and Eq. 1. Val- 
ues of uj were taken from the X-ray Data Booklet 16 . For 
Ni and Ca, /kc* = 0-89 and 0.88, respectively, assuming 
neutral atoms (Bambynek et al. 1972). We assume /kq 
= 0.9 for all other elements. The values of EW ca \ c , calcu- 
lated assuming solar abundances and a covering fraction 
/c=1.0, are listed in Table 3. 

Relative abundances are calculated as, e.g., Zsi/Zp e 
= (jEW r ob s f Si/-BW r ca i C| si) / (J5W r obs l Fe/-EWcai c> Fe), where 
E Wcaic,Fe has been estimated using the abundance derived 
from the Fe K edge depth, Zp e =l.l7, and the solar abun- 
dance was used for £W ca i C) si- We caution that for Ar, Ca, 
and Ni, the lines are weak and the statistics are poor. For 
Si, there may be systematic effects due to current XIS cal- 
ibration as well as from the fact that the continuum near 
1.7 keV is mixture of all three power-law components. Ta- 
ble 5 lists the abundance ratios. There is agreement on 
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the values of Zs / Zp e derived from the edges and from the 
lines. 

We observe [m/H] = log(Zpe) = +0.1, a slightly higher 
metallicity value than in the metal-rich population of glob- 
ular clusters at 2-20 kpc radii ([m/H] ~ —0.1, Peng et al. 
2004) or in the outer stellar bulge component (e.g., [m/H] 
~ -0.2 at a radius of 8 kpc, Harris & Harris 2002). It 
is only slightly higher than values for the ISM of other 
early- type galaxies (Humphrey Sz Buote 2006). Our ob- 
served value of Zug/Z-pe is consistent with the value of 2 
obtained by Peng et al. (2004) for Cen A’s globular clus- 
ters. It is similar to values for the stellar components of 
large ellipticals in general, [m/H] ~ 0.0 — 0.4 (Henry Sz 
Worthey 1999). Trager et al. (2000) also found the stellar 
populations of 8 field ellipticals to have [m/H] ~ 0.0 — 0.4, 
though it may not be straightforward to compare Cen A 
to field ellipticals due to Cen A’s recent merger. 

However, our observed [m/H] value is much higher than 
the average value found in the halo stars of Cen A, [m/H] 
= -0.4 at radii of 20-30 kpc (Harris & Harris 2000, Harris, 
Harris Sz Poole 1999). It is therefore unlikely that the cir- 
cumnuclear material and the diffuse plasma have directly 
originated in the relatively metal-poor outer halo, unless 
enrichment via local star formation has occurred. 

The pre-merger origin of the circumnuclear material 
cannot be known for certain, though it is plausible that 
during the merger, gaseous material originally in the outer 
portions of the progenitor galaxies lost sufficient angu- 
lar momentum to be transported to the nuclear regions 
(Toomre Sz Toomre 1972). Furthermore, connections be- 
tween the nuclear gas and starburst activity at radii of 
hundreds of pc and less are well-known (e.g., Hopkins et 
al. 2006 and references therein). The likelihood that Cen 
A’s circumnuclear gas has been enriched due to starburst 
processes is therefore high. 

The degree of metallicity enrichment depends on many 
factors, including the IMF, star formation rate, and the 
total gas consumed, but simulations of disk-disk mergers 
incorporating star formation by Mihos Sz Hernquist (1994) 
demonstrated that [m/H] values of 0.3-0. 4 at the center of 
the resulting merged ellipticals are plausible. The metal- 
licity value derived from the Suzaku data is consistent with 
enrichment due to the ongoing star formation in the inner 
bulge, as a result of the merging process. 

The observed relative element abundances also support 
enrichment. We calculated the expected relative abun- 
dances due to enrichment by Type II and Type la SNe, 
ignoring enrichment via stellar wind processes. We used 
the abundances relative to solar as found in Tsujimoto et 
al. (1995; their Figures 1-2), assuming a Salpeter IMF, 
converting to the abundances of Lodders (2003), and as- 
suming one Type la explosion for every ten Type II explo- 
sions. For the ratios listed in Table 5, we expect Zo/Zp e 
~2.5, Zyjsje / Zp e through Z$\/Zp e to be 1.4-1. 7, and Z$/Zp e 
through Zca/^Fe and Z^i/Zp e to all be ~ 0.7-0. 9. Most of 
our observed abundance ratios are consistent (or at least 
roughly consistent) with these predictions. The current 
burst of star formation is estimated to have started at 
least 50 Myr ago (e.g., Dufour et al. 1979). Assuming 
star formation at a uniform rate of 30 M© yr -1 (Telesco 
1978 measured a star formation rate ten times that of the 
Milky Way’s spiral arms), this is sufficient to enrich the 


gas to the current metallicity levels within that time span 
(M. Loewenstein, priv. comm.), although we cannot know 
for certain how much enrichment of the gas has occurred 
within the last 10 8 yr and how much occurred closer in 
time to the merger. 

5. CONCLUSIONS 

The combination of Suzaku ' s XIS and HXD has allowed 
us to obtain a high quality spectrum of the nucleus (inner 
2 kpc) of Cen A spanning 0.3-250 keV. The long expo- 
sure time, narrow response of the XIS, and the low back- 
ground make the XIS spectrum one of the highest qual- 
ity CCD spectra of Cen A to date, particularly below 2 
keV. The HXD-PIN above 12 keV is more sensitive than 
RXTE-BEXTE or INTEGRAL instruments, allowing us 
to accurately and directly measure the 12-76 keV flux and 
spectral shape. The HXD-GSO provides a direct estimate 
of the 45-250 keV flux of Cen A. 

Our best-fit model includes several broadband compo- 
nents. The >3 keV emission can be fit by two absorbed 
power-laws. Including data down to 0.3 keV, we find that 
two VAPEC components plus a third absorbed power-law 
component can fit the data. The soft X-rays show emis- 
sion lines due to He- and H-like O, Fe XVII, He- and 
H-like Ne, and He-like Mg. Two VAPEC components of 
temperatures k&T = 0.2 and 0.6 keV can explain the soft 
X-ray emission lines. The third power-law component rep- 
resents a blend of diffuse emission and knot emission from 
the kpc-scale jet, plus emission from a population of XRBs 
and other X-ray point sources, plus a likely contribution 
from scattered nuclear emission, which is likely less than 
1% of the total hard X-ray power-law emission. 

In the hard X-rays, the primary power-law is absorbed 
by a column of about 1.5 x 10 23 cm 2 . The secondary power- 
law, with a brightness about 19% of the primary, is ab- 
sorbed by a much higher column, about 7xl0 23 cm 2 . In 
our best-fit model, the photon indices were tied, and the 
best fit value was 1.817to!oio- Untying the photon in- 
dices did not result in a significant fit improvement. The 
Suzaku data are consistent with two scenarios previously 
forwarded to explain the nature of the dual power-laws. 
As suggested by Evans et al. (2004), one power-law may 
be due to emission from the sub-pc VLBI jet, while the 
other may be associated with Bondi accretion at the black 
hole, though it is not obvious which corresponds to the pri- 
mary/secondary power-law. Alternatively, a partial cover- 
ing model, as suggested by Turner et al. (1997), may apply: 
84% (16%) of the sky as seen from the X-ray source is ob- 
scured by a column 1.5 (7) xlO 23 cm 2 . In any event, it 
more likely that the material obscuring the hard X-rays is 
associated with the circumnuclear material than with the 
optical dust lane. 

The hard X-ray photon index, the 2-10 keV flux, the 
Fe K a line intensity and the column density of the ma- 
terial absorbing the primary power-law are all consistent 
with recent observations over the last ~5 years, suggesting 
that no large variations in these parameters have occurred 
over the last few years. During the Suzaku observation, 
the hard X-ray flux was not strongly variable on short 
timescales, showing only a 10% increase in flux through- 
out the observation. 

K-shell absorption edges due to Fe, Ca, and S are signif- 
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icantly detected. Several fluorescent emission lines, likely 
from the absorbing material in transmission, are detected. 
In addition to Fe K a and Fe K/3, Ka lines from Si, S, 
Ar, Ca and Ni are detected at at least 90% confidence. 
The latter three are detected for the first time in Cen A. 
The strict upper limit on the amount of Compton reflec- 
tion, R < 0.05, supports the notion that the fluorescent 
lines originate in Compton-thin material. The Fe Ka line 
width is < 2500 km s -1 FWHM. If the line-emitting mate- 
rial is in virial orbit, then it is no closer than 200 light-days 
to the black hole. No emission due to He- or H-like Fe is 
detected. The Fe K edge energy is consistent with ab- 
sorption by Fe I; ionization stages above Fe ~V are ruled 
out. The gas does not likely completely surround the X- 
ray continuum source; the equivalent width of the Fe Ka 
line is more consistent with gas which is clumpy and/or 
covers the continuum source only partially, or the bulk of 
the gas may lie off the line of sight to the X-ray source, 
possibly in the form of a disk. 

Metallicities are measured for the circumnuclear absorb- 
ing gas for the first time in Cen A. The high signal- to-noise 
ratio of the XIS spectrum has enabled us to use the K-shell 
absorption edge depths, relative fluorescent line strengths, 
and thermal plasma emission lines to derive abundances 
(we have assumed abundances in the thermal plasma gas 
are same as those for the circumnuclear absorbing mate- 
rial). We observe a value of [m/H] = log(Zpe) = +0.1 for 


the circumnuclear material/diffuse plasma which is much 
higher than for the halo stars at 20-30 kpc radii. It is likely 
that the merger which created Cen A’s famous edge-on 
dust disk also triggered the infall of gas to the inner re- 
gions of Cen A, where it currently accretes onto the black 
hole. Our metallicity observation suggests that the ac- 
creting material could not have originated in the relatively 
metal-poor outer halo unless enrichment due to local star 
formation occurred at some point. The relative observed 
abundances are indeed consistent with enrichment by a 
mixture of Type II and Type la supernovae. 

This high signal-to-noise CCD observation is just a taste 
of what can be done with an orbiting X-ray calorimeter. In 
particular, Const ellation-X will be able to provide abun- 
dance constraints on the circumnuclear material in poten- 
tially dozens of AGNs once it is launched. 

The authors gratefully acknowledge the dedication and 
hard work of the Suzaku hardware teams and operations 
staff for making this observation possible and for assis- 
tance with data calibration and analysis. A.M. thanks 
Rick Rothschild, Mike Loewenstein and Tim Kallman 
for useful discussions. This research has made use of 
HEASARC online services, supported by NASA/GSFC. 
This research has also made use of the NASA/IPAC Extra- 
galactic Database, operated by JPL/California Institute of 
Technology, under contract with NASA. 
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Table 1 

Soft X-ray Emission Lines 


Line 

Observed Energy 

Intensity 

EW 

Identification 

Centroid (keV) 

(10 -5 ph cm -2 s -1 ) 

(eV) 

O VII 

n ££>1+0.006 
0.564_o 0 04 

10.9 ± 2.1 

39±0.8 

O VIII 

0.654±0.003 

9*3 + i'i 

37±\° 

Fe XVII 3s-2p / O VII RRC 

n 7*30+0.004 

u. /oz_0 006 

£ 3+1- 0 

O.O_0 7 

23l i 

Fe XVII 3d-2p ( x Pi) / Fe XVII 3d-2p ( 3 £>i) 

0.820±g;gg| 


33^6 

Ne IX 

0.901±0.004 

6.3tg;J 

36+| 

NeX 

1.010±0.003 

4.1±0.6 

33±5 

Mg XI 

■f oQi +0.008 

l.ooi_Q 004 

1.5±0.4 

22±6 


Note. — Equivalent widths EW (Col. [4]) were determined relative to a locally-fit continuum. 


Table 2 


Best-fit parameters for selected models 


Parameter 

Model 5 

Model 6 

Model 7 

Model 8 

X 2 / do f 

9806.8/9X58 

9745.1/9159 

9639.6/9144 

9529.9/9142 

r i j2 

1 qco+0.008 

±.ooz_o 019 

1 QQ 1 +0.019 
i.ooi _0 oio 

i oon+0-016 

i.ozu_o oio 

1 Qi 7+0.023 
I.oi. i — 0.010 

PL 1 Norm . 1 

0.116iO°? 2 

0 . 112 ± 0.002 

0 . 111 db 0.002 

0 111 + 0-001 
0.002 

Nh,i (/ 10 22 cm 2 ) 

i5-2iSi 

15.1±0.2 

14.7±0.2 

14 7+0 - 3 
i4 - — 0.2 

PL 2 Norm . 1 

0 09fi+0’0°4 

u.uzo_o 002 

0.022T0.003 

0 . 02 li °;°° 2 

0.021±0.003 

Ah , 2 (/ 10 22 cm 2 ) 

82±7 

75±8 

72l 10 

70+ 11 

VAPEC 1 k B T (keV) 

n qi +0.03 
u - ,5 - L — 0.02 

0.24±0.02 

0.23±0.02 

0 . 22 l °;° 2 

VAPEC 1 Norm. 

7.2^i><10 - 4 

5.0l^xl0 - 4 

SOioiJxlO -4 

4.3l°; 1 xio - 4 

VAPEC2 k B T (keV) 

(1 74+009 
u< ‘ 0.03 

0.62±0.01 

0.62±0.01 

0.56i«;°3 

VAPEC2 Norm. 

3.7^;| xlO " 4 

5.1±0.3xl0 -4 

5.2l^|xl0- 4 

4.8±0.4xl0 ~ 4 

r 3 

1 16+0- 13 
1 * i 0 - 0 .01 

1-28 toil 

1 41+ 014 

1 qi+0-08 

PL3 Norm . 1 

6.9 tol x 10~ 4 

7-4 toi x 10 - 4 

7.8l^| x 10- 4 

7.01q3 x 10-4 

Nh ,3 (/ 10 22 cm 2 ) 

0 (fixed) 

<0.03 

<0.05 

<0.03 


Note. — Best-fit parameters for selected broadband models. All models include two absorbed X-ray power-laws (PL1, PL2), 
two VAPEC components, a third soft X-ray power-law (PL3), and Fe Kq and K (3 emission lines. Listed are: Model 5, where 
PL3 is not absorbed; Model 6, where PL3 is absorbed; Model 7, where ZVPHABS abundances are thawed for S, Ca, and Fe and 
fluorescent lines for Si, S, Ar, Ca and Ni are added; and Model 8, our best-fit model, in which the VAPEC abundances Zfq e = 
2.7dt0.3 and ^Mg = 1.9i0.4. 

1 : units are ph keV -1 cm -2 s _1 at 1 keV. 


Table 3 


Fluorescent Emission Line Parameters 


Line 

Energy 

Centroid (keV) 

Intensity (10 5 
ph cm -2 s -1 ) 

EW ohs 

(eV) 

Ax 2 

F-test 

Prob. 

EW cal c 
(eV) 

Si Kq 

1.71±0.01 

0.9±0.2 

24±6 

39.8 

7.6x 10 _1 ° 

26.3 

S Kq 

2.307±0.016 

0.8±0.4 

6±2 

5.0 

9.6x 10 -2 

11.4 

Ar Kq 

2.994±0.023 

1.2±0.5 

13±4 

20.1 

8.0x10-* 

4.7 

Ca Kq 

3.690±0.023 

1.8±0.7 

8±3 

6.2 

5.4xl0- 2 

3.7 

Fe Kq 

6.394±0.003 

23.2±1.0 

83±3 

3077 

<lxlO- 40 

128 

FeK (3 

7.056 1 

3.0 2 

ll 2 

28.3 

2.7x 10 -7 

17 2 

Ni Kq 

7.47±0.05 

1.4±0.8 

15±5 

9.3 

l.lxlO -2 

7.1 


Note. — Results are for Model 8. Observed equivalent widths (Col. [4]) are determined relative to a locally- fit continuum. 
All widths were tied to that for the Fe Kq line, 14+ eV. EW ca \ c (Col. [7]) denotes predicted equivalent width values using 
Eq. 1 (uniform covering assuming line of sight column density), solar abundances, and a covering fraction of 1.0. 

1 denotes a fixed parameter. 

2 denotes value tied to 0.13 that for the Fe Kq line. 
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0 5x1 0 4 10 5 

Time (sec since 2005 Aug 19, 04:11 UT) 

Fig. 1.- — Orbitally binned XIS and HXD light curves. The top panel shows the 2-10 keV count rate light curve, summed over all four XIS 
cameras. The 12-76 keV HXD-PIN light curve is shown in the middle panel. The 45-250 keV HXD-GSO light curve is shown in the bottom 
panel. 




10 20 50 50 100 200 

Observed Energy (keV) Observed Energy (keV) 

Fig. 2. — HXD-PIN (left) and GSO (right) spectra. The upper panels show the net source spectrum (gray points ), the background (lower 
black points ), and the total (source + background) spectrum (upper black points). The PIN spectra have been binned such that the net 
spectrum has a minimum signal-to-noise ratio of 10a per bin. The GSO spectra have been binned as described in §2. The lower panels shows 
the ratio of the net source spectrum to the total spectrum. 
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Table 4 


Hard X-ray Absorption Edges 


K-shell 

Edge 

ZVPHABS 

Abundance 

Edge Energy 
(keV) 

Edge Optical 
Depth, r 

S 

Ca 

Fe 

1 14+ 0 - 11 

i c+0.8 
i *°-0.6 

i ' i '-0.09 

2.468+g;°n 

4.040zb0.045 

7.103±0.015 

0 24 - * -0 ' 03 
U-^_0.06 

0.03 zb 0.01 

0.21 ± 0.01 


Note. — The ZVPHABS abundances in Col. (2) are from Model 8. The edge energy and optical depths (cols. [3]- [4]) are from 
a model with ZVPHABS abundances for Fe, Ca and S set to zero and edges fit instead. 


Table 5 


Abundance Ratios 


Method 

Ratio 

Value 

1 

Zo/Zpe 

0.8 ± 0.4 

1 

Z Ne /Z Fe 

2.3 zb 0.3 

1 

Zug/Z^e 

1.6 ± 0.3 

3 

Zsi/ Zj? e 

1.5 zb 0.4 

2 

Zs / Zp e 

1.0 ± 0.1 

3 

Zs/Z-pe 

0.8 ± 0.3 

3 

Z A.T / Zp Q 

4 ± 1 

2 

Zca/^Fe 

1.3 dz 0.6 

3 

ZCa/^Fe 

3 ± 1 

3 

Zm/Zp G 

3 zb 1 


Note. — Methods 1, 2, and 3 (Col. [1]) denote abundances derived from thermal plasma emission lines, K shell edge depths, 
and Ka fluorescent emission line intensities, respectively. 



10 100 
Energy (keV) 

Fig. 3. — Residuals to Model 1, the single absorbed power-law (top, and Model 2, which included dual absorbed power-law components and 
Fe Ka and K/3 line emission (bottom). Black points denote the three FI XIS spectra, which have been co-added here for clarity. Gray, red 
and blue points denote XIS-1, PIN, and GSO, respectively. Rest- frame energies are shown. The XIS data have been plotted with a binning 
factor of 10. 




14 


Markowitz et al. 



0.5 1 1.5 

Energy (keV) 

Fig. 4. — Ratio of the soft band data to a simple power-law, showing the prominent emission lines. Black points denote the XIS BI; gray 
points denote data from XIS-0, 2 and 3, which have been co-added for clarity. Rest-frame energies are shown. The data have been plotted 
with a binning factor of 3. 
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Energy (keV) 

Fig. 5. — Residuals to various broadband models are shown. From top to bottom are: Model 3, with soft X-rays fit with a simple power-law; 
Model 4, with soft X-rays fit with two VAPEC components only; Model 6, with soft X-rays fit with two vapec components plus an absorbed 
power-law; Model 7, with K-shell absorption edges and Ka emission lines for S, Si, Ar, Ca and Ni fit; and Model 8, our best-fit model, with 
VAPEC abundances for Ne and Mg thawed. Black points denote the three FI XIS spectra, which have been co-added here for clarity. Gray, 
red and blue points denote XIS-1, PIN, and GSO, respectively. Rest-frame energies are shown. The XIS data have been plotted with a 
binning factor of 10. 



r PIN/XIS2 norm. 

Fig. 6. — Contour plots of the Compton reflection fraction R versus the hard X-ray photon index ri=r 2 and the PIN/XIS-2 instrument 
normalization for Model 8. Dotted, dashed, and solid lines denote 68, 95.4, and 99.73% confidence levels, respectively. 
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6 7 8 

Energy (keV) 

Fig. 7. — Residuals to a simple power-law model in the Fe K bandpass, showing the prominent Fe Ka emission line as well as the Fe K/3 
line at 7.08=1=0.03 keV, the Ni Ka line at 7.47d=0.05 keV, and the Fe K edge. Black points denote the three FI XIS spectra, which have been 
co-added for clarity; gray points denote XIS-1. Rest- frame energies are shown. All data have been plotted with a binning factor of 10. 



Energy (keV) 

Fig. 8. — Residuals to a simple power-law model in the S K and Si K bandpass, showing the Si Ka emission line at 1.71=1=0.01 keV, the 
S Ka line at 2.307=1=0.016 keV, and the S K edge. Black points denote the three FI XIS spectra, which have been co-added for clarity; gray 
points denote XIS-1. Rest- frame energies are shown. All data have been plotted with a binning factor of 10. 



3 4 5 

Energy (keV) 

Fig. 9. — Residuals to a simple power-law model in the Ar K and Ca K bandpass, showing the Ar Ka emission line at 2.994d=0.023 keV, 
the Ca Ka line at 3.690±0.023 keV, and the Ca K edge. Black points denote the three FI XIS spectra, which have been co-added for clarity; 
gray points denote XIS-1. Rest-frame energies are shown. All data have been plotted with a binning factor of 10. 
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1.75 1.8 1.85 1.9 0.09 0.1 0.11 0.12 

T PL Norm. 

(ph keV -1 cm -2 s -1 at 1 keV) 


Fig. 10. — Contour plots of the line of sight column density absorbing the primary power-law versus the photon index (left) and 1 keV 
normalization (right) of the primary power-law component in Model 8. Dotted, dashed, and solid lines denote 68, 95.4, and 99.73% confidence 
levels, respectively. The plus signs mark the best-fit values. 



Observed Energy (keV) 

Fig. 11. — Unfolded broadband spectrum for Model 8, illustrating the three absorbed power-law components and the two vapec components. 
All data have been plotted with a binning factor of 10. 





